SELP Asp603Asn and severe thrombosis in COVID-19 males.

Thromboembolism is a frequent cause of severity and mortality in COVID-19. However, the etiology of this phenomenon is not well understood. A cohort of 1186 subjects, from the GEN-COVID consortium, infected by SARS-CoV-2 with different severity was stratified by sex and adjusted by age. Then, common coding variants from whole exome sequencing were mined by LASSO logistic regression. The homozygosity of the cell adhesion molecule P-selectin gene (SELP) rs6127 (c.1807G > A; p.Asp603Asn) which has been already associated with thrombotic risk is found to be associated with severity in the male subcohort of 513 subjects (odds ratio = 2.27, 95% Confidence Interval 1.54-3.36). As the SELP gene is downregulated by testosterone, the odd ratio is increased in males older than 50 (OR 2.42, 95% CI 1.53-3.82). Asn/Asn homozygotes have increased D-dimers values especially when associated with poly Q ≥ 23 in the androgen receptor (OR 3.26, 95% CI 1.41-7.52). These results provide a rationale for the repurposing of antibodies against P-selectin as adjuvant therapy in rs6127 male homozygotes especially if older than 50 or with an impaired androgen receptor.

Gain- and Loss-of-Function CFTR Alleles Are Associated with COVID-19 Clinical Outcomes.

Carriers of single pathogenic variants of the CFTR (cystic fibrosis transmembrane conductance regulator) gene have a higher risk of severe COVID-19 and 14-day death. The machine learning post-Mendelian model pinpointed CFTR as a bidirectional modulator of COVID-19 outcomes. Here, we demonstrate that the rare complex allele [G576V;R668C] is associated with a milder disease via a gain-of-function mechanism. Conversely, CFTR ultra-rare alleles with reduced function are associated with disease severity either alone (dominant disorder) or with another hypomorphic allele in the second chromosome (recessive disorder) with a global residual CFTR activity between 50 to 91%. Furthermore, we characterized novel CFTR complex alleles, including [A238V;F508del], [R74W;D1270N;V201M], [I1027T;F508del], [I506V;D1168G], and simple alleles, including R347C, F1052V, Y625N, I328V, K68E, A309D, A252T, G542*, V562I, R1066H, I506V, I807M, which lead to a reduced CFTR function and thus, to more severe COVID-19. In conclusion, CFTR genetic analysis is an important tool in identifying patients at risk of severe COVID-19.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative.

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

An explainable model of host genetic interactions linked to COVID-19 severity.

We employed a multifaceted computational strategy to identify the genetic factors contributing to increased risk of severe COVID-19 infection from a Whole Exome Sequencing (WES) dataset of a cohort of 2000 Italian patients. We coupled a stratified k-fold screening, to rank variants more associated with severity, with the training of multiple supervised classifiers, to predict severity based on screened features. Feature importance analysis from tree-based models allowed us to identify 16 variants with the highest support which, together with age and gender covariates, were found to be most predictive of COVID-19 severity. When tested on a follow-up cohort, our ensemble of models predicted severity with high accuracy (ACC = 81.88%; AUCROC = 96%; MCC = 61.55%). Our model recapitulated a vast literature of emerging molecular mechanisms and genetic factors linked to COVID-19 response and extends previous landmark Genome-Wide Association Studies (GWAS). It revealed a network of interplaying genetic signatures converging on established immune system and inflammatory processes linked to viral infection response. It also identified additional processes cross-talking with immune pathways, such as GPCR signaling, which might offer additional opportunities for therapeutic intervention and patient stratification. Publicly available PheWAS datasets revealed that several variants were significantly associated with phenotypic traits such as "Respiratory or thoracic disease", supporting their link with COVID-19 severity outcome.

Multiomic analysis reveals cell-type-specific molecular determinants of COVID-19 severity.

The determinants of severe COVID-19 in healthy adults are poorly understood, which limits the opportunity for early intervention. We present a multiomic analysis using machine learning to characterize the genomic basis of COVID-19 severity. We use single-cell multiome profiling of human lungs to link genetic signals to cell-type-specific functions. We discover >1,000 risk genes across 19 cell types, which account for 77% of the SNP-based heritability for severe disease. Genetic risk is particularly focused within natural killer (NK) cells and T cells, placing the dysfunction of these cells upstream of severe disease. Mendelian randomization and single-cell profiling of human NK cells support the role of NK cells and further localize genetic risk to CD56bright NK cells, which are key cytokine producers during the innate immune response. Rare variant analysis confirms the enrichment of severe-disease-associated genetic variation within NK-cell risk genes. Our study provides insights into the pathogenesis of severe COVID-19 with potential therapeutic targets.

Pathogen-sugar interactions revealed by universal saturation transfer analysis.

Many pathogens exploit host cell-surface glycans. However, precise analyses of glycan ligands binding with heavily modified pathogen proteins can be confounded by overlapping sugar signals and/or compounded with known experimental constraints. Universal saturation transfer analysis (uSTA) builds on existing nuclear magnetic resonance spectroscopy to provide an automated workflow for quantitating protein-ligand interactions. uSTA reveals that early-pandemic, B-origin-lineage severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike trimer binds sialoside sugars in an "end-on" manner. uSTA-guided modeling and a high-resolution cryo-electron microscopy structure implicate the spike N-terminal domain (NTD) and confirm end-on binding. This finding rationalizes the effect of NTD mutations that abolish sugar binding in SARS-CoV-2 variants of concern. Together with genetic variance analyses in early pandemic patient cohorts, this binding implicates a sialylated polylactosamine motif found on tetraantennary N-linked glycoproteins deep in the human lung as potentially relevant to virulence and/or zoonosis.

Carriers of ADAMTS13 Rare Variants Are at High Risk of Life-Threatening COVID-19.

Thrombosis of small and large vessels is reported as a key player in COVID-19 severity. However, host genetic determinants of this susceptibility are still unclear. Congenital Thrombotic Thrombocytopenic Purpura is a severe autosomal recessive disorder characterized by uncleaved ultra-large vWF and thrombotic microangiopathy, frequently triggered by infections. Carriers are reported to be asymptomatic. Exome analysis of about 3000 SARS-CoV-2 infected subjects of different severities, belonging to the GEN-COVID cohort, revealed the specific role of vWF cleaving enzyme ADAMTS13 (A disintegrin-like and metalloprotease with thrombospondin type 1 motif, 13). We report here that ultra-rare variants in a heterozygous state lead to a rare form of COVID-19 characterized by hyper-inflammation signs, which segregates in families as an autosomal dominant disorder conditioned by SARS-CoV-2 infection, sex, and age. This has clinical relevance due to the availability of drugs such as Caplacizumab, which inhibits vWF-platelet interaction, and Crizanlizumab, which, by inhibiting P-selectin binding to its ligands, prevents leukocyte recruitment and platelet aggregation at the site of vascular damage.

Host genetic basis of COVID-19: from methodologies to genes.

The COVID-19 pandemic caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is having a massive impact on public health, societies, and economies worldwide. Despite the ongoing vaccination program, treating COVID-19 remains a high priority; thus, a better understanding of the disease is urgently needed. Initially, susceptibility was associated with age, sex, and other prior existing comorbidities. However, as these conditions alone could not explain the highly variable clinical manifestations of SARS-CoV-2 infection, the attention was shifted toward the identification of the genetic basis of COVID-19. Thanks to international collaborations like The COVID-19 Host Genetics Initiative, it became possible the elucidation of numerous genetic markers that are not only likely to help in explaining the varied clinical outcomes of COVID-19 patients but can also guide the development of novel diagnostics and therapeutics. Within this framework, this review delineates GWAS and Burden test as traditional methodologies employed so far for the discovery of the human genetic basis of COVID-19, with particular attention to recently emerged predictive models such as the post-Mendelian model. A summary table with the main genome-wide significant genomic loci is provided. Besides, various common and rare variants identified in genes like TLR7, CFTR, ACE2, TMPRSS2, TLR3, and SELP are further described in detail to illustrate their association with disease severity.

The polymorphism L412F in TLR3 inhibits autophagy and is a marker of severe COVID-19 in males.

The polymorphism L412F in TLR3 has been associated with several infectious diseases. However, the mechanism underlying this association is still unexplored. Here, we show that the L412F polymorphism in TLR3 is a marker of severity in COVID-19. This association increases in the sub-cohort of males. Impaired macroautophagy/autophagy and reduced TNF/TNFα production was demonstrated in HEK293 cells transfected with TLR3L412F-encoding plasmid and stimulated with specific agonist poly(I:C). A statistically significant reduced survival at 28 days was shown in L412F COVID-19 patients treated with the autophagy-inhibitor hydroxychloroquine (p = 0.038). An increased frequency of autoimmune disorders such as co-morbidity was found in L412F COVID-19 males with specific class II HLA haplotypes prone to autoantigen presentation. Our analyses indicate that L412F polymorphism makes males at risk of severe COVID-19 and provides a rationale for reinterpreting clinical trials considering autophagy pathways.Abbreviations: AP: autophagosome; AUC: area under the curve; BafA1: bafilomycin A1; COVID-19: coronavirus disease-2019; HCQ: hydroxychloroquine; RAP: rapamycin; ROC: receiver operating characteristic; SARS-CoV-2: severe acute respiratory syndrome coronavirus 2; TLR: toll like receptor; TNF/TNF-α: tumor necrosis factor.

Rare variants in Toll-like receptor 7 results in functional impairment and downregulation of cytokine-mediated signaling in COVID-19 patients.

Toll-like receptors (TLR) are crucial components in the initiation of innate immune responses to a variety of pathogens, triggering the production of pro-inflammatory cytokines and type I and II interferons, which are responsible for innate antiviral responses. Among the different TLRs, TLR7 recognizes several single-stranded RNA viruses including SARS-CoV-2. We and others identified rare loss-of-function variants in X-chromosomal TLR7 in young men with severe COVID-19 and with no prior history of major chronic diseases, that were associated with impaired TLR7 signaling as well as type I and II IFN responses. Here, we performed RNA sequencing to investigate transcriptome variations following imiquimod stimulation of peripheral blood mononuclear cells isolated from patients carrying previously identified hypomorphic, hypofunctional, and loss-of-function TLR7 variants. Our investigation revealed a profound impairment of the TLR7 pathway in patients carrying loss-of-function variants. Of note, a failure in IFNγ upregulation following stimulation was also observed in cells harboring the hypofunctional and hypomorphic variants. We also identified new TLR7 variants in severely affected male patients for which a functional characterization of the TLR7 pathway was performed demonstrating a decrease in mRNA levels in the IFNα, IFNγ, RSAD2, ACOD1, IFIT2, and CXCL10 genes.

Common, low-frequency, rare, and ultra-rare coding variants contribute to COVID-19 severity.

The combined impact of common and rare exonic variants in COVID-19 host genetics is currently insufficiently understood. Here, common and rare variants from whole-exome sequencing data of about 4000 SARS-CoV-2-positive individuals were used to define an interpretable machine-learning model for predicting COVID-19 severity. First, variants were converted into separate sets of Boolean features, depending on the absence or the presence of variants in each gene. An ensemble of LASSO logistic regression models was used to identify the most informative Boolean features with respect to the genetic bases of severity. The Boolean features selected by these logistic models were combined into an Integrated PolyGenic Score that offers a synthetic and interpretable index for describing the contribution of host genetics in COVID-19 severity, as demonstrated through testing in several independent cohorts. Selected features belong to ultra-rare, rare, low-frequency, and common variants, including those in linkage disequilibrium with known GWAS loci. Noteworthily, around one quarter of the selected genes are sex-specific. Pathway analysis of the selected genes associated with COVID-19 severity reflected the multi-organ nature of the disease. The proposed model might provide useful information for developing diagnostics and therapeutics, while also being able to guide bedside disease management.

Severe COVID-19 in Hospitalized Carriers of Single CFTR Pathogenic Variants.

The clinical presentation of COVID-19 is extremely heterogeneous, ranging from asymptomatic to severely ill patients. Thus, host genetic factors may be involved in determining disease presentation and progression. Given that carriers of single cystic fibrosis (CF)-causing variants of the CFTR gene-CF-carriers-are more susceptible to respiratory tract infections, our aim was to determine their likelihood of undergoing severe COVID-19. We implemented a cohort study of 874 individuals diagnosed with COVID-19, during the first pandemic wave in Italy. Whole exome sequencing was performed and validated CF-causing variants were identified. Forty subjects (16 females and 24 males) were found to be CF-carriers. Among mechanically ventilated patients, CF-carriers were more represented (8.7%) and they were significantly (p < 0.05) younger (mean age 51 years) compared to noncarriers (mean age 61.42 years). Furthermore, in the whole cohort, the age of male CF-carriers was lower, compared to noncarriers (p < 0.05). CF-carriers had a relative risk of presenting an abnormal inflammatory response (CRP ≥ 20 mg/dL) of 1.69 (p < 0.05) and their hazard ratio of death at day 14 was 3.10 (p < 0.05) in a multivariate regression model, adjusted for age, sex and comorbidities. In conclusion, CF-carriers are more susceptible to the severe form of COVID-19, showing also higher risk of 14-day death.

Association of Toll-like receptor 7 variants with life-threatening COVID-19 disease in males: findings from a nested case-control study.

Background: Recently, loss-of-function variants in TLR7 were identified in two families in which COVID-19 segregates like an X-linked recessive disorder environmentally conditioned by SARS-CoV-2. We investigated whether the two families represent the tip of the iceberg of a subset of COVID-19 male patients. Methods: This is a nested case-control study in which we compared male participants with extreme phenotype selected from the Italian GEN-COVID cohort of SARS-CoV-2-infected participants (<60 y, 79 severe cases versus 77 control cases). We applied the LASSO Logistic Regression analysis, considering only rare variants on young male subsets with extreme phenotype, picking up TLR7 as the most important susceptibility gene. Results: Overall, we found TLR7 deleterious variants in 2.1% of severely affected males and in none of the asymptomatic participants. The functional gene expression profile analysis demonstrated a reduction in TLR7-related gene expression in patients compared with controls demonstrating an impairment in type I and II IFN responses. Conclusions: Young males with TLR7 loss-of-function variants and severe COVID-19 represent a subset of male patients contributing to disease susceptibility in up to 2% of severe COVID-19. Funding: Funded by private donors for the Host Genetics Research Project, the Intesa San Paolo for 2020 charity fund, and the Host Genetics Initiative. Clinical trial number: NCT04549831.

Shorter androgen receptor polyQ alleles protect against life-threatening COVID-19 disease in European males.

BACKGROUND: While SARS-CoV-2 similarly infects men and women, COVID-19 outcome is less favorable in men. Variability in COVID-19 severity may be explained by differences in the host genome. METHODS: We compared poly-amino acids variability from WES data in severely affected COVID-19 patients versus SARS-CoV-2 PCR-positive oligo-asymptomatic subjects. FINDINGS: Shorter polyQ alleles (≤22) in the androgen receptor (AR) conferred protection against severe outcome in COVID-19 in the first tested cohort (both males and females) of 638 Italian subjects. The association between long polyQ alleles (≥23) and severe clinical outcome (p = 0.024) was also validated in an independent cohort of Spanish men <60 years of age (p = 0.014). Testosterone was higher in subjects with AR long-polyQ, possibly indicating receptor resistance (p = 0.042 Mann-Whitney U test). Inappropriately low serum testosterone level among carriers of the long-polyQ alleles (p = 0.0004 Mann-Whitney U test) predicted the need for intensive care in COVID-19 infected men. In agreement with the known anti-inflammatory action of testosterone, patients with long-polyQ and age ≥60 years had increased levels of CRP (p = 0.018, not accounting for multiple testing). INTERPRETATION: We identify the first genetic polymorphism that appears to predispose some men to develop more severe disease. Failure of the endocrine feedback to overcome AR signaling defects by increasing testosterone levels during the infection leads to the polyQ tract becoming dominant to serum testosterone levels for the clinical outcome. These results may contribute to designing reliable clinical and public health measures and provide a rationale to test testosterone as adjuvant therapy in men with COVID-19 expressing long AR polyQ repeats. FUNDING: MIUR project "Dipartimenti di Eccellenza 2018-2020" to Department of Medical Biotechnologies University of Siena, Italy (Italian D.L. n.18 March 17, 2020) and "Bando Ricerca COVID-19 Toscana" project to Azienda Ospedaliero-Universitaria Senese. Private donors for COVID-19 research and charity funds from Intesa San Paolo.

Employing a systematic approach to biobanking and analyzing clinical and genetic data for advancing COVID-19 research.

Within the GEN-COVID Multicenter Study, biospecimens from more than 1000 SARS-CoV-2 positive individuals have thus far been collected in the GEN-COVID Biobank (GCB). Sample types include whole blood, plasma, serum, leukocytes, and DNA. The GCB links samples to detailed clinical data available in the GEN-COVID Patient Registry (GCPR). It includes hospitalized patients (74.25%), broken down into intubated, treated by CPAP-biPAP, treated with O2 supplementation, and without respiratory support (9.5%, 18.4%, 31.55% and 14.8, respectively); and non-hospitalized subjects (25.75%), either pauci- or asymptomatic. More than 150 clinical patient-level data fields have been collected and binarized for further statistics according to the organs/systems primarily affected by COVID-19: heart, liver, pancreas, kidney, chemosensors, innate or adaptive immunity, and clotting system. Hierarchical clustering analysis identified five main clinical categories: (1) severe multisystemic failure with either thromboembolic or pancreatic variant; (2) cytokine storm type, either severe with liver involvement or moderate; (3) moderate heart type, either with or without liver damage; (4) moderate multisystemic involvement, either with or without liver damage; (5) mild, either with or without hyposmia. GCB and GCPR are further linked to the GCGDR, which includes data from whole-exome sequencing and high-density SNP genotyping. The data are available for sharing through the Network for Italian Genomes, found within the COVID-19 dedicated section. The study objective is to systematize this comprehensive data collection and begin identifying multi-organ involvement in COVID-19, defining genetic parameters for infection susceptibility within the population, and mapping genetically COVID-19 severity and clinical complexity among patients.

Clinical and molecular characterization of COVID-19 hospitalized patients.

Clinical and molecular characterization by Whole Exome Sequencing (WES) is reported in 35 COVID-19 patients attending the University Hospital in Siena, Italy, from April 7 to May 7, 2020. Eighty percent of patients required respiratory assistance, half of them being on mechanical ventilation. Fiftyone percent had hepatic involvement and hyposmia was ascertained in 3 patients. Searching for common genes by collapsing methods against 150 WES of controls of the Italian population failed to give straightforward statistically significant results with the exception of two genes. This result is not unexpected since we are facing the most challenging common disorder triggered by environmental factors with a strong underlying heritability (50%). The lesson learned from Autism-Spectrum-Disorders prompted us to re-analyse the cohort treating each patient as an independent case, following a Mendelian-like model. We identified for each patient an average of 2.5 pathogenic mutations involved in virus infection susceptibility and pinpointing to one or more rare disorder(s). To our knowledge, this is the first report on WES and COVID-19. Our results suggest a combined model for COVID-19 susceptibility with a number of common susceptibility genes which represent the favorite background in which additional host private mutations may determine disease progression.

ACE2 gene variants may underlie interindividual variability and susceptibility to COVID-19 in the Italian population.

In December 2019, an initial cluster of interstitial bilateral pneumonia emerged in Wuhan, China. A human-to-human transmission was assumed and a previously unrecognized entity, termed coronavirus disease-19 (COVID-19) due to a novel coronavirus (SARS-CoV-2) was described. The infection has rapidly spread out all over the world and Italy has been the first European country experiencing the endemic wave with unexpected clinical severity in comparison with Asian countries. It has been shown that SARS-CoV-2 utilizes angiotensin converting enzyme 2 (ACE2) as host receptor and host proteases for cell surface binding and internalization. Thus, a predisposing genetic background can give reason for interindividual disease susceptibility and/or severity. Taking advantage of the Network of Italian Genomes (NIG), here we mined whole-exome sequencing data of 6930 Italian control individuals from five different centers looking for ACE2 variants. A number of variants with a potential impact on protein stability were identified. Among these, three more common missense changes, p.(Asn720Asp), p.(Lys26Arg), and p.(Gly211Arg) were predicted to interfere with protein structure and stabilization. Rare variants likely interfering with the internalization process, namely p.(Leu351Val) and p.(Pro389His), predicted to interfere with SARS-CoV-2 spike protein binding, were also observed. Comparison of ACE2 WES data between a cohort of 131 patients and 258 controls allowed identifying a statistically significant (P value < 0.029) higher allelic variability in controls compared with patients. These findings suggest that a predisposing genetic background may contribute to the observed interindividual clinical variability associated with COVID-19, allowing an evidence-based risk assessment leading to personalized preventive measures and therapeutic options.